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的 重要 内 容 与 基础 。 早 期 定量 生物 测定 数据 的 整理 
分 析 多 采用 Finney 机 值 分 析 方 法 , 但 因 其 时 间 与 剂 
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摘要 : 生物 测定 是 生物 学 、 医 学 、 毒 理学 的 重要 内 容 与 基础 。 常 用 的 定量 生物 测定 数据 分 析 方 法 时 间 - 剂 量 -死亡 率 
模型 (TDM) 不 能 对 复杂 生 测 数据 建立 统一 模型 , 信息 利用 不 充分 。 本 文 基于 文 持 癌 量 回归 (SVR), 提出 了 一 种 能 对 
不 同 供 试 因子 、 不 同 供 试 对 象 和 不 同 环境 条 件 下 复杂 生 测 数据 统一 建 模 的 新 方法 。14 个 简单 生 测 数据 和 2 套 复杂 
生 测 数据 的 对 比分 析 绪 果 表 明 ,SVR 模型 拟 合 与 留 一 法 预测 精度 均 优 于 TDM 模型 , 估计 的 LDso。 和 LT; 等 指标 更 为 
可 信 。SVR 模型 有 望 作为 TDM 模型 的 有 益 补 充 , 在 定量 生物 测定 数据 分 析 中 得 到 广泛 应 用 。 
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Abstract: Bioassay plays an important role in the studies of biology, medicine and toxicology. The time- 
dose-mortality model (TDM ) widely applied to quantitative bioassay data analysis can not construct a 
unified model for complex bioassay data, and has the disadvantage of utilizing the information incompjetely. 
Based on support vector regression (SVR), a novel quantitative bioassay model has been developed which 
can construct a unified model for complex data with different test factors, different test objects and different 
environment factors. We compared the prediction performance between SVR and TDM using 14 simple data 
and 2 complex data. The results showed that SVR achieved better precision than TDM not only in self- 
consistency test but also in jackknife test, implying that the estimated values of LD;, and LT;, by the former 
are more reliable. As a useful supplement to TDM, SVR has the potential to be widely used for quantitative 
bioassay data analysis. 
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Robertson ，1989 ;， Robertson and Preisler，1992 ) ， 并 
由 冯 明 光 等 (1996 ) 介绍 引 进 国 内 ， 它 综合 考虑 了 时 
间 与 剂量 效应 并 将 其 统一 到 同一 个 模型 中 , 对 多 数 


量 效应 的 评价 相互 脱节 、 不 能 揭示 时 间 与 剂量 的 互 
作 效 应 、 未 充分 利用 实验 信息 等 缺陷 , 近年 来 已 逐 
步 被 时 间 - 剂 量 -死亡 率 模 型 (time-dose-mortality 
model, TDM ) 所 取代 。TDM 按 其 数学 结构 又 称 为 互 
补 重 对 数 模型 ( complementary log-log model，CLL ) ， 
最 早 由 Preisler 和 Robertson 提出 (Preisler and 


简单 生物 测定 数据 具有 较 好 的 拟 合 能 力 , 并 方便 给 
出 LD; 和 LTs 等 各 种 指标 ,目前 已 广泛 应 用 。 
然而 , TDM 模型 仍然 存在 以 下 葬 端 : (1) CILL 
是 一 个 具 显 性 表达 式 的 非 线性 经 验 模型 ,很 难 假定 
所 有 生物 测定 数据 均 符 合同 一 确定 性 模型 而 仅 参 数 
变化 。 对 CLL 模型 拟 合 优 度 的 评估 常 采 用 Pearson 
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卡 方 检验 或 Hosmer-Lemeshow 检验 ( 唐 局 义 和 汉 明 
光 , 2007) 。 在 实际 应 用 中 ,即使 采用 较 不 严格 的 
Hosmer-Lemeshow 检验 , 仍 有 部 分 生 测 数据 不 能 通 
过 。 排 除 实验 误差 造成 的 影响 , 对 这 一 事实 的 最 可 
能 解释 是 : 这 些 生 测 数据 本 吴 不 适 于 用 CLL 模型 分 
析 。(2) 昆虫 死亡 与 否 属 二 项 分 布 , 对 CLL 模型 参 
数 的 鸽 计 在 DPS 软件 中 是 通过 二 项 分 布 最 大 似 然 
画 数 逼 近来 获得 的 , 但 实际 应 用 中 有 些 生 测 数 据 不 
符合 二 项 分 布 (如 杀菌 剂 或 除草 剂 的 时 间 - 剂 量 - 相 
对 抑制 率 )， 限制 了 CLL 模型 的 广泛 应 用 ( 柏 连 阳 
等 , 2003) 。 最 近 , 唐 尘 等 (2006 ) 采用 非 线 性 最 小 
二 乘法 估计 CLL 模型 参数 , 提出 了 杀菌 剂 生 物 测定 
动态 评价 时 间 - 剂 量 - 抑 菌 率 模型 , 但 该 参数 估计 方 
法 及 由 此 获得 的 LD;, 和 LT; 等 指标 估计 精度 尚 需 
更 多 研究 实例 的 支持 。(3) 基于 CLL 的 时 间 - 剂 量 - 
死亡 率 模 型 与 时 间 - 剂 量 - 抑 菌 率 模 型 不 能 对 不 同 供 
试 因子 、 不 同 供 试 对 象 和 不 同 环境 下 的 复杂 生物 测 
定数 据 建立 统一 模型 ， 导 致 实验 信息 利用 不 充分 。 
例如 分 析 在 m 种 温度 下 nn 个 球 孢 白 僵 亲 Beauveria 
bassiana 菌株 不 同 剂 量 对 桃 蚜 Myzus persicae 不 同时 
间 段 的 死亡 率 , 需 构建 m x n 个 独立 的 CLL 模型 。 
我 们 知道 , 统计 分 析 中 多 样本 的 比较 需 采 用 方差 分 
析 , 如 对 多 样本 采用 上 上 检验 进行 两 两 比较 , 不 但 繁 
琐 、 犯 a 错误 的 可 能 性 增 大 , 估计 试验 误差 时 的 精 
度 也 有 所 损失 。 类 似 地 , 对 一 组 复杂 生物 测定 数据 
也 应 统一 到 同一 个 模型 中 , 而 不 宜 构 建 多 个 独立 的 
CLL 模型 。 为 方便 表述 , 本文 称 单一 供 试 因 子 、 单 
一 供 试 对 象 和 单一 环境 下 的 生物 测定 数据 为 简单 数 
据 , 多 种 供 试 因 子 和 /或 多 种 供 试 对 象 和 /或 多 种 环 
境 下 的 生物 测定 数据 为 复杂 数据 。(4) TDM 模型 
隶属 于 Fisher 经 典 参 数 统计 理论 , 把 回归 分 析 表 达 
为 特定 参数 化 模型 (CLL 模型 ) 的 参数 估计 问题 ; 但 
经 典 参数 统计 理论 是 样本 数目 趋 于 无 穷 大 时 的 渐 近 
理论 , 它 基 于 经 验 风 险 最 小 (训练 误差 最 小 ) , 要求 
训练 样本 充分 大 。 实 际 生物 测定 数据 往往 是 小 样本 
(这 里 样本 数 指 的 是 剂量 与 时 间 段 等 形成 的 组 合 
数 , 而 不 是 每 一 剂量 下 的 供 试 生物 个 体 数 )。 在 有 
限 样 本 下 ， 训 练 误差 小 并 不 总 能 获得 好 的 预测 效 
采 , 某 些 情况 下 , 训练 误差 过 小 反而 会 导致 推广 能 
力 的 下 降 , 即 真 实 风 险 的 增加 , 这 就 是 过 拟 合 或 过 
学 习 问 题 (Vapnik , 1995)。(5) 生物 测定 数据 分 析 
的 主要 目标 仅仅 是 求 出 未 知 函 数 在 目标 点 的 值 ， 如 
LD 和 LTs。 经 典 参 数 统计 理论 采用 归纳 和 演绎 分 
两 步 解 决 : 第 一 步 (从 特殊 到 一 般 , 即 归 纳 ), 利用 


最 大 似 然 法 从 一 个 给 定 的 函数 集 (CLL 模型 ) 解 决 
一 个 一 般 性 问题 ; 第 二 步 ( 从 一 般 到 特殊 ， 即 演 
绎 ) ,从 一 般 性 问题 出 发 , 计算 目标 点 的 值 。 当 拥 
有 信息 量 ( 小 样本 ) 只 够 解决 欲求 问题 , 却 不 足以 解 
决 一 个 一 般 性 问题 时 ,用 一 个 给 定 的 函数 集 来 耳 接 
佑 计 目 标点 上 的 函数 值 , 形成 了 一 种 新 的 从 特殊 到 
特殊 的 推理 方法 一 一 转 导 推理 (Vapnik, 1995 ) 。 转 
导 推 理 即 小 样本 集 推理 的 主要 原则 是 : 如 果 对 欲求 
解 的 某 一 问题 只 拥有 有 限 的 信息 , 那么 应 该 直接 求 
解 问题 , 而 决 不 能 求解 一 个 更 一 般 的 问题 来 作为 一 
个 中 间 步 又， 因为 可 能 所 拥有 的 信息 足以 直接 求解 
问题 , 但 不 足以 解决 一 个 更 一 般 的 中 间 问 题 ( 邓 旋 
扬 和 田 英 杰 , 2004 ) 。 

支持 癌 量 机 (support vector machine ,SVM ) 是 统 
计 学 习 理 论 的 集大成 者 , 它 基 于 结构 风险 最 小 与 转 
导 推 理 , 较 好 地 解决 了 小 样本 、 非 线性 、 过 拟 合 、 
维 数 灾难 和 局 部 极 小 等 问题 , 泛 化 推广 能 力 优异 
(Vapnik，1995 ; 邓 力 扬 和 田 英 杰 , 2004) 。SVM 包 
括 支 持 回 量 分 类 (support vector classification, SVC) 
和 支持 回 量 回归 (support vector regression，SVR ) 。 
本 文 针 对 TDM 模型 存在 的 缺陷 , 建立 了 基于 SVR 
的 定量 生物 测定 数据 分 析 新 方法 , 并 以 多 个 实例 验 
证 了 新 方法 的 有 效 性 。 


1 数据 与 方法 


1.1 数据 来 源 

14 个 简单 生物 测定 数据 和 2 个 复杂 生物 测定 
数据 引 目 表 1 和 2 所 列 文献 。 
1.2 LIBSVM 2.8 软件 包 与 SVR 

LIBSVM 2.8 软件 包 的 SVR 功能 需 用 到 4 个 常 
用 程序 : svmscale 用 于 对 原始 数据 规格 化 ，svmtrain 
用 于 训练 ，svmpredict 用 于 预测 ，gridregression. py 
用 于 自动 搜索 核 呜 数 最 优 参 数 c, g,p (ce[-1， 
6],gel[l-8,0],pe[-8, -1 ，, 步 长 均 为 1)。 
各 程序 用 法 及 其 参数 设置 参见 Chang 和 Lin 
(2001 ) 。 
1.3 基于 SVR 的 生 测 数据 分 析 
1.3.1 拟 合 : 以 观察 累积 死亡 率 y(% ) 为 因 变 量 ， 
供 试 因子 x 、 供 试 对 象 x,、 剂 量 x，、 时 间 x 、 环 
境 因 子 ( 如 温度 ) xs 等 为 目 变 量 构成 二 维 数据 表 。 
其 中 供 试 因子 x 与 供 试 对 象 %, 为 单位 回 量 。 如 供 试 
对 和 象 为 小 菜 蛾 Plutella xylostella 对 Bt 敏感 品系 与 抗 
性 品系 , 则 敏感 品系 (4 %,1 ) 可 编码 为 (0, 1) ， 抗 性 
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有 后 系 ( 2,2 ) 编码 为 (1,;0) 。 一 般 地 , 复杂 生 测 数据 
中 供 试 因子 和 供 试 对 和 象 回 量 维 数 分 别 为 供 试 因 子 数 
和 供 试 对 象 类 别 数 。 人 简单 生 测 数据 目 变 量 不 包括 供 
试 因子 、 供 斌 对象 与 环境 因子 , 仅 含 剂量 与 时 间 。 

数据 集 经 svmscale 规格 化 后 , 选用 径 回 基 核 函 
数 , 以 gridregression. py 经 格 点 搜索 获得 最 优 参数 
c, 8,P 了 组合 。 取 最 优 参 数组 合 ， 以 svmtrain 训练 建 
模 ， 以 svmpredict 预测 数据 集 目 喇 ,， 获得 拟 合 均 方 
误差 (root mean square error，RMSE ) 与 决定 系数 
( determination coefficient, R*') ,作为 模型 拟 合 优 劣 
的 评价 标准 。 





(1) 


ZL7(i,true) 一 Jane | y(i,predict) 一 Yediet ] 





ES[y(i,tue) -ye] SLyCi, predict) 一 yaedic]” 
(2) 

式 中 , y(i, true) 为 第 i 个 累积 死亡 率 观察 值 ， 
y(i,，predict) 为 第 i 个 累积 死亡 率 预 测 值 , y,, 为 真 
实 值 均值 , yweas 为 预测 值 均值 。 
1.3.2 留 一 法 预测 : 模型 评估 不 仅 希 望 拟 合 精度 
高 , 更 重要 的 是 泛 化 推广 能 力 强 。 可 采用 留 一 法 
(Leave-One-Out method ，L00 ) 检验 SVR 模型 的 独 
立 预 测 能 力 : 从 nn 个 样本 中 每 次 拿 出 一 个 样本 作为 
测试 集 , 以 剩 下 的 -1 个 样本 构成 训练 集 , 依次 经 
gridregression. py，svmtrain 和 svmpredict 等 完成 预 
测 。 同 上 计算 RMSE 与 R 值 。 
1.3.3 估计 LD;, 和 LTs: 在 给 定 供 试 因 子 、 供 试 
对 象 、 环 境 因 了 于 和 时 间 下 求 LD;, 可 产生 一 个 仅 剂 
量 按 小 步 长 变化 的 测试 集 , 以 1.3.1 训练 好 的 SVR 
模型 预测 ,预测 累积 死亡 率 最 接近 0.5 所 对 应 的 剂 
量 即 LDx。 类 似 可 求 LTso。 
1.4 参 比 模型 TDM 

TDM 模型 的 累积 死亡 率 拟 合 值 、LD; 与 LTsx 由 
DPS 软件 直接 给 出 ( 唐 启 义 和 冯 明光 , 2007) ,并 据 
此 算得 拟 合 RMSE 与 R 值 。DPS 软件 不 允许 原始 
数据 单个 缺失 , 因此 基于 DPS 不 能 直接 给 出 TDM 
模型 留 一 法 预测 结果 。 根 据 TDM 模型 原理 
( Robertson and Preisler, 1992 ) , 本 文 以 Matlab 语言 
重新 编号 了 TDM 留 一 法 预测 程序 并 在 Linux 上 运 
行 通过 。 经 非 缺 失 数 据 运算 检验 ,结果 与 DPS 一 
致 。 对 一 套 复 杂 生 测 数据 , 先 以 多 个 独立 的 CLL 和 模 


型 分 别 实施 留 一 法 预测 , 最 后 合并 预测 结 
留 一 法 预测 RMSE 与 R 值 。 


2 结果 与 分 析 


基于 SVR 的 14 个 简单 生 测 数 据 和 2 个 复杂 生 
测 数据 定量 分 析 结 果 见 表 1 和 2。 从 RMSE 值 与 尺 
值 看 , 无 论 是 简单 生 测 数据 还 是 复杂 生 测 数据 ,也 
无 论 是 模型 拟 合 还 是 留 一 法 预测 ，SVR 模型 均 优 于 
TDM 模型 。 


, 计算 


3 讨论 


定量 生物 测定 数据 分 析 的 目标 是 估计 LDx 和 
LT 等 指标 以 评估 供 试 因子 对 供 试 对 象 的 毒 力 大 
小 , 其 实质 是 独立 预测 。 在 小 样本 条 件 下 , 传统 基 
于 经 验 风险 最 小 的 回归 模型 的 主要 缺陷 是 往往 模型 
拟 合 精度 尚 可 , 但 独立 预测 精度 明显 下 降 。 因 此 ， 
评估 模型 更 应 注重 其 独立 预测 能 力 。 本 文 涉 及 的 
14 个 简单 生 测 数 据 和 2 套 复杂 生 测 数据 是 检索 有 
完整 原始 数据 文献 后 随机 选取 的 , 结果 显示 SVR 
模型 无 论 是 拟 合 还 是 留 一 法 预测 精度 均 优 于 TDM 
模型 ， 因 而 有 理由 相信 SVR 模型 估计 的 LDx 和 
LT 等 指标 更 为 可 靠 。 

TDM 模型 对 不 同 供 试 因 子 、 不 同 供 斌 对象 和 不 
同 环境 下 的 复杂 生物 测定 数据 需 单独 建 模 , 因而 实 
验 信 息 利 用 不 充分 。 例 如 ,， 某 球 孢 日 僵 菌 对 桃 虹 
在 20C 和 25C 下 的 生 测 数据 分 析 , 假定 观测 时 间 
为 处 理 后 1 -10 d 共 10 个 时 间 单 位 、 球 孢 日 僵 坎 
对 数 浓度 为 1 ~7 共 7 个 浓度 , 假设 20 时 第 10 
天 的 LCa 徘 近 对 数 浓度 5，TDM 建 模 估计 207Y 的 
LCso 时 采用 20 时 第 1 天 与 对 数 浓度 1 的 观察 数 
据 , 而 不 采用 25 必 时 第 10 天 与 对 数 浓度 5 的 观察 
数据 , 直觉 上 就 欠 合 理 。 相 反 ，SVR 模型 可 对 不 
同 供 试 因子 、 不 同 供 试 对 象 和 不 同 环境 条 件 下 复 
杂 生 测 数据 统一 建 模 ,其 实验 信息 利用 更 为 充分 ， 
且 可 估计 未 观测 环境 条 件 下 (如 22 ) 的 LDs、 
LTs 等 指标 。 

SVR 定量 生物 测定 分 析 模 型 的 主要 缺陷 包括 : 
(1) SVR 是 非 参数 模型 , 而 TDM 是 机 理 模型 , 后 者 
模型 中 各 个 参数 有 明确 的 生物 学 意义 。(2) SVR 
没有 考虑 试验 样本 大 小 , 也 不 能 给 出 LD 和 LT; 等 
估计 指标 的 置信 区 间 。 相 反 , TDM 建 模 时 将 每 头 试 
忠 作为 一 个 样本 进行 了 考虑 ， 日 能 给 出 LD;。 和 LT 
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等 佑 计 指 标的 置信 区 间 。(3) SVR 应 用 稍 嫌 复 杂 ， ” 现 阶 段 , SVR 模型 可 作为 定量 生物 测定 数据 分 析 中 
运算 速度 偏 慢 。 而 TDM 使 用 简便 , 运算 速度 快 。 TDM 模型 的 一 个 有 益 补 充 。 


表 1 简单 生 测 数 据 的 SVR 模型 和 TDM 模型 比较 
Table 1 Comparison of SVR model and IDM model on simple bioassay data 


SVR 模型 TDM 模型 
供 试 因子 ” 供 试 对 象 SVR model TDM moqel 参考 文献 
Factors Objects 拟 合 Fitting 留 一 法 L00 拟 合 Fitting 留 一 法 LO0 References 
LD LT LD. LT 
RMSE (R*) RMSEL(R’) RMSE (R*) RMSEL(R) 
0. 0263 0.0412 0. 0693 0.0849 Robertson and 
Ac Oh 0.00 1.9 -0.51 1.5 
(0.9904) (0.9761) (0.9330) (0. 8978) Preisler, 1992 
0. 0447 0.0557 0. 0520 0.0574 刘 彩 玲 和 冯 明 
Za Mp 1.58 4.9 1.53 4.6 7 
(0.9702) (0.9536) (0.9594 (0.9504) 光 , 1998 
» y 0.0094 0.0140 jy 65 38 0. 0249 0.0271 0.60 3 地 和 民生 2003 
n 。 。 。 。 当 ， 
(0.9988) (0.9974) (0.9916) (0.9900) Ms 
0. 0922 0.0964 0.0964 0. 1034 鲍 佳 生 和 冯 明 
Za Mp 1.47 4.6 1.60 4.5 人 
(0.8768) (0.8671) (0. 8593) (0. 8385) 光 , 2006 
0. 0520 0. 0583 0.0557 0. 0663 许 谦和 冯 明 
Zr Mp 1.48 2.7 1.48 2.6 + 说 
(0.9706) (0.9612) (0.9653) (0.9504) 光 , 1999 
pd 0.0279 0.0312 0. 0332 0. 0346 
Nl 1.28 6.9 1.46 7.5 Xu etal., 1999 
(F95127) (0.9675) (0.9592) (0.9528) (0.9479) 
pd 0.0138 0.0187 0. 0490 0.0510 
Nl 1.38 3.5 1.41 5.0 Xu eal., 1999 
(F95129) (0.9962) (0.9932) (0.9543) (0.9485) 
0.0277 0.0374 0. 0447 0. 0490 
Pd Mp 0.32 3.2 0.29 3.1 梁 勇 , 2003 
(0.9906) (0.9829) (0.9752) (0.9708) 
0.0195 0. 0346 0. 0387 0. 0458 
Ct Mp 1.84 4.0 1.80 4.5 李 伟 , 2002 
(0.9928) (0.9777) (0.9712) (0.9604) 
Bb 0.0181 0. 0294 0. 0458 0. 0500 
Tt 5.30 4.8 5.01 4.7 徐 华 潮 ,2003 
(F263) (0.9964) (0.9900) (0.9759) (0.9708) 
Bb 0.0361 0. 0458 0. 0500 0. 0566 
Tt 4.70 4.1 4.89 3.9 徐 华 潮 ，2003 
( 浙 B) (0.9864) (0.9785) (0.9738) (0.9667) 
Bb 0.0193 0.0332 0.0361 0.0412 
Tt 4.30 3.9 3.99 3.7 徐 华 潮 ，2003 
( B12) (0.9970) (0.9912) (0.9890) (0.9860) 
Bb 0. 0447 0. 0707 0.0806 0. 0883 
Tt 4.40 3.9 4.52 4.2 徐 华 潮 ,2003 
( REBbO1) (0.9799) (0.9495) (0.9339) (0.9203) 
Bb 0. 0230 0. 0436 0. 0447 0. 0529 
Tt 4.80 3.6 4.38 3.5 徐 华 潮 ，2003 
( Bxs) (0.9956) (0.9843) (0.9833) (0.9773) 


LDso : 仅 给 出 最 长 观察 时 间 间 隔 下 的 LDso 或 LCso Only the LDso or LCso according to the maximum time interval was listed; LTs : 仅 给 出 最 高 剂量 
或 最 高 浓度 下 的 LTso Only the LTso according to the maximum dose or the highest concentration was listed; Ac: 乙酰 甲 胺 磷 Acephate; Za: 安徽 虫 辣 
霉 Zoophthora anhuiensis; Pn: 诺 氏 虫 疙 霉 Pandora nouryi; Zr: 根 虫 瘟 霉 Zoophthora radicans; Pd: 飞 乱 虫 六 霉 Pandora delphacis; Ct: 块 耳 霉 
Conidiobolus thromboides ; Bb: 球 隐 白 僵 菌 Beauveria bassiana; Oh: 0.，herbivorous; Mp: 桃 蚜 Myzus persicae; Nl: 褐飞虱 Nilaparvata lugens; Tt: 竹 
梢 凸 层 斑 蚜 Takecallis taiwanus， 下 辣 The same below. 
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表 2 复杂 生 测 数据 的 SVR 模型 和 TDM 模型 比较 
Table 2 Comparison of SVR model and TDM model on complex bioassay data 


SVR 模型 
供 试 因子 供 试 对 象 VR moe 
Factors Objects 拟 合 Fitting 留 一 法 LOO 
LD50 LTso 
RMSE(R2) RMSE(R) 
Pd 
1.38 7.3 
(F95127) NI 0.0197 0.0346 
pd (0.9908) (0.9710) 
1.48 3.5 
(F95129) 
Bb 
5.30 4.7 
(上 -263 ) 
Bb 
4.70 4.1 
( 浙 B) 
Bb 0.0346 0.0490 
Tt 4.00 3.9 
( B12) (0.9890) (0.9661) 
Bb 
4.20 3.8 
( REBbO01 ) 
Bb 
4.50 3.5 
( Bxs) 
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